Universality and the three-body parameter of helium-4 trimers 
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We consider a system of three helium-4 atoms, which is so far the simplest realistic three-body 
system exhibiting the Efimov effect, in order to analyse deviations from the universal Efimov three- 
body spectrum. We first calculate the bound states using a realistic two-body potential, and then 
analyse how they can be reproduced by simple effective models beyond Efimov's universal theory. 
We find that the non-universal variations of the first two states can be well reproduced by models 
parametrized with only three quantities: the scattering length and effective range of the original 
potential, and the strength of a small three-body force. Furthermore, the three-body parameter 
which fixes the origin of the infinite set of three-body levels is found to be consistent with recent 
experimental observations in other atomic species. 



I. INTRODUCTION 

The universal attraction found by V. Efimov [T] 
for any quantum system of three particles interacting 
through short-range interactions with a large scatter- 
ing length has now been evidenced in many exper- 
iments using ultra-cold atoms pHl7j. In particular, 
Efimov trimers, i.e. three-body bound states resulting 
from this attraction, were observed as a function of the 
scattering length. Because these trimers are unnat- 
urally large compared with the range of the interac- 
tions, they have universal properties determined solely 
by the universal attraction and a few parameters. In 
particular, their spectrum has a simple structure with 
discrete scale invariance. This was confirmed experi- 
mentally to some degree, but there appeared quanti- 
tative deviations from this structure. 

One of the main reasons is the fact that at most the 
first two states of the spectrum could be observed so 
far. It is known that these first states do not follow 
accurately the universal behaviour expected for the 
higher excited states (more loosely bound states), be- 
cause their size is not very large compared to the range 
of the interactions and therefore they still depend on 
the details of these interactions. However, it is quite 
involved to correctly describe these interactions for 
three atoms because of their complex hyperfine struc- 
ture and the lack of knowledge of the three-body po- 
tential surfaces. For this reason, experimental results 
have been interpreted so far using either interaction- 
dependent corrections to the universal theory of Efi- 
mov [Pol 116], fT8| [19] , or by other effective models re- 
producing the two-body physics in the energy range 
of the observed trimers |20If2"2~] . 

While these effective approaches could reproduce 
the experimental results to some extent, it remains un- 
clear on a theoretical basis why they could do so. For 
example, corrections to the Efimov universal theory 
sometimes required to introduce an ad-hoc variation 
of a 3-body parameter to explain the data |15[ [TBI [22] . 
Another puzzling fact is that effective two-body model 
approaches could reproduce some of the deviations 
from universal theory observed in the experimental 
data, suggesting that they could be explained by two- 
body interactions only [2DH2"!2] . While it is established 



that in general the knowledge of two-body interac- 
tions only is not enough to accurately determine the 
energy of Efimov trimers |23| . the contribution from 
realistic two-body interactions to the short-range 3- 
body phase and non-universal deviations are not fully 
understood. 

The purpose of this paper is to clarify these issues 
by testing the effective approaches with respect to 
the numerically exact solution of a realistic theoretical 
model. We choose 4 He3, as it is the simplest triatomic 
system with van der Waals inteactions which exhibits 
the Efimov attraction jHHH]. 

The paper is organised as follows. In section II, we 
review some of the effective models used to describe 
Efimov trimer experiments. In section III, we present 
realistic calculations for 4 He 3 , and how they are re- 
produced by these effective models. 



II. EFFECTIVE MODELS FOR EFIMOV 
PHYSICS 

A. Efimov's universal theory 

The essence of the Efimov effect is the appearance 
of an effective attractive potential —sfjR 2 attrac- 
tion between three particles with very large scattering 
length. Here, so is a number approximately equal to 
1.00624 for identical bosonic particles, and R denotes 
the hyper-radius (average distance between particles) , 

R^l^+rh+rl), (1) 

where 7"i2,r23, and r%i are the three particles' relative 
distances. This attraction can lead to the existence of 
three-body bound states, the so-called Efimov trimers. 
Because it is a long-range attraction, trimers with 
sufficiently small binding energy extend to large dis- 
tances where the interactions are negligible. The only 
effect of interactions is to set short-distance boundary 
conditions on the three-body wave function. The first 
boundary condition occurs when two particles come 
close to one another, but within a distance r larger 
than the range of their interaction. There, the wave 
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function ip has to satisfy the Bethe-Peierls boundary 
condition 
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where a is the s-wave scattering length of the two- 
body interaction, which fixes the phase of the two- 
body wave function accumulated from short distance 
at low energy. The second boundary condition occurs 
when three particles come close together, but still at 
distances R larger than the range of their interactions. 
The wave function has to satisfy the Efimov boundary 
condition 
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where A is the so-called Efimov three-body parameter, 
which fixes the accumulated phase of the three-body 
wave function at low energy. 

The Efimov theory thus relies on only two parame- 
ters, a and A, to fully describe the three-body physics 
in the low-energy and large a regime. When nor- 
malised in units of A, the trimer energy spectrum ex- 
hibits a universal structure as a function of aA, rep- 
resented in Fig. [5] There is an accumulation point 
in the spectrum at a = oo, that is when the two- 
body interaction is resonant, and the whole spectrum 
is invariant by a discrete scale transformation, namely 
multiplying all distances by e n ^ s ° ~ 22.7, which fol- 
lows from Eq. (|3|. The wave functions and energies 
can be calculated numerically by either solving the 
3-body free Schrodinger equation with conditions Q 
and ([3| , or equivalently solving its corresponding inte- 
gral equation which is known as the Skorniakov- Ter- 
Martirosian equation |29| . 
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where F is the unknown function (related to the full 
wave function *£), and e = t# is the renormalised 
energy E. Here, the upper bound P of the integral 
sets the three-body phase, and can be related to the 
Efimov three-body parameter by |22| : 

. , arctan sq + im . 

P = Aexp( ) 

so 



B. Non- universal models 



where k = s/mE/h. The universal limit corresponds 
to the first term, which is set by the zero-energy scat- 
tering length a. The next order defines the effective 
range r e . It is straightforward to generalise Eq. Q 
to the non-universal regime by replacing the zero- 
energy scattering length a by the energy-dependent 
scattering length a(E) |22l [30] . Likewise, the cutoff 
P is expected to be replaced by an energy-dependant 
quantity P(E) [55]. With these replacements, Eq. Q 
corresponds to the most general contact model with 
energy-dependent boundary conditions. 

Although the energy dependence of a(E) is gener- 
ally known, that of P(E) is presently unknown. It is 
one of the purposes of this paper to investigate this 
dependence from comparison with other models. 



2. Two-body effective interaction 

Another approach is to replace the real interaction 
by a simple effective interaction with the same low- 
energy spectrum. One possible choice is a Gaussian 
potential [3T1I32"] : 
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which is parametrised by Vq and tq to reproduce both 
the scattering length and effective range. Another 
convenient choice is a separable interaction |20ff22] : 



v = v \<t>)(4>\ 



(7) 



where the state \<f>) can also be chosen to be a Gaussian 
function </>(r) = e^^ r ^' ^ for simplicity. The advantage 
of separable potentials is that they have formal simi- 
larities with contact interactions, and lead to a simple 
integral equation similar to Eq. (HI: 
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~ p2) F{q) = 



where G(x) = exp[J^° ^-j-dt]. One can check that the 
integrand of (JsJ) tends to that of (|4j) at low momenta, 
but decays at high momenta q > l/»"o, which removes 
the need to introduce an upper bound to the integral. 



1. Non-universal corrections 



III. REALISTIC AND EFFECTIVE 
CALCULATIONS FOR 4 He 3 



The deviations of two-body physics from universal- 
ity at low energy are well known and can be encoded in 
the energy variation of the two-body phase shift 6(E), 
or equivalently a two-body scattering length a(E). At 
low energy, we have the following effective-range low- 
energy expansion: 
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A. Realistic calculations with LM2M2 potentials 

To model 4 He interactions realistically, we choose 
the LM2M2 potential [33] to describe the two-body 
interactions. This potential has a repulsive hard core 
at short distance, and a van der Waals tail —Ce/r G 
at large distance, as shown in Fig. [T] Its scatter- 
ing length is 100.01 A which is 18.6 times the van 



3 



10 



5 - 




2 4 6 8 

Interatomic distance (A) 

Figure 1: LM2M2 potential |33| used for the realistic cal- 
culations. The dotted curve indicates the van der Waals 
asymptote —C%jr®. 



der Waals length l v( \W = (mCe/h) 1 ^ 4 . The three- 
body interaction has been shown to bring only small 
corrections |34fl5rj] . and we neglect it in this study. 
Thus, in our calculation, the 3-body phase which fixes 
the energy of Efimov trimers builds up only from the 
LM2M2 two-body interaction. 

The 3-body Schrodinger equation with the LM2M2 
potential is solved numerically using the Gaussian ex- 
pansion method (GEM). This method was proposed 
as a means to perform accurate calculation for three- 
and four-body systems (37] . In this method, a well- 
chosen set of Gaussian basis function is used, forming 
an approximately complete set in a finite coordinate 
space, so that one can describe accurately both short- 
range correlation and the long-range asymptotic be- 
haviour of the wavefunction for bound systems as well 
as for scattering states. It was demonstrated that the 
GEM provides the same caliber of numerical precision 
as, for example, the Faddeev-Yakubovsky method for 
3 H ( 3 He) and 4 He, and can be used to address various 
kinds of few-body problems in atomic, baryonic and 
quark-level systems [37]. 

In order to solve three-body 4 He trimer problem, 
we use three sets of Jacobi coordinates illustrated in 
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Figure 2: Jacobi coordinates for all the rearrangement 
channels (c = 1 ~ 3) of 4 He trimer system. Three 4 He 
atoms are to be symmetrized. 

The Schrodinger equation and the total Hamilto- 
nian are given by 

(H - E) % JM = , (9) 



H = T + V(n) + V(r 2 ) + V(r 3 ) , (10) 



where T is the kinetic-energy operator and V(ri), 
V{r2) and V{r$) are the interactions between two 4 He 
atoms, described by the LM2M2 potential V(r). 

The total three-body wavefunction is described as 
a sum of amplitudes for all rearrangement channels 
(c = 1 ~ 3) of Fig. § 

e=l n,N t,L 

(11) 

We take the functional forms of <fi n i m (r), i/inlm(R-) 

as 

<Wr) = r l e-<r' r "?Y lm (r) , 

iPnlm (R) = R L e-W R ^Y LM (R) , (12) 

where the Gaussian range parameters are chosen ac- 
cording to geometrical progressions: 

r n = ria™" 1 (n = 1, . . . , n max ) , 
Rn = R\A N ~ l (N=l,...,N maK ) . (13) 

The eigenenergies E in Eq|9]and the coefficients C in 
Eq|TT]are determined by the Rayleigh-Ritz variational 
method. 

Although the scattering length of 4 He is already 
large compared to the range of its two-body poten- 
tial V(r), we did several calculations for rescaled 
potentials XV (r) in order to cause a divergence of 
the scattering length, as was done in previous stud- 
ies |241 1231 12"7] . This enables us to mimic the broad 
Feshbach resonances used in ultra-cold atom experi- 
ments [35] , and better appreciate the Efimov structure 
of the spectrum. The scattering length a diverges for 
A = 0.97412, while the physical results for real 4 He 
correspond to A = 1. As A is varied near the di- 
vergence of a, the effective range r e also changes but 
always remains positive and on the order of the scaled 
van der Waals length A 1/,4 Z vc [-yy, and the ratio r e /a re- 
mains a monotonic function of A, as shown in Fig. [3J 
For this reason, we choose to report our results as a 
function of the scattering length in units of the effec- 
tive range, rather than A itself. 

The LM2M2 potential supports one two-body 
bound state and its energy variation with the scat- 
tering length is represented in Fig. [4] When the scat- 
tering length becomes smaller than its physical value 
(A > 1.0), this dimer energy significantly deviates 
from the universal limit of small binding energy and 
large scattering length (Eq. ( |l4"| ) with a(E) —> a). The 
results for three atoms are shown in Fig. [5] Our re- 
sults are in very good agreement with the most accu- 
rate calculations using the LM2M2 potential [391 140|. 
One can see that the first two trimers' energies quali- 
tatively follow Efimov's universal spectrum, but as in 
the two-body case, there are significant deviations for 
small scattering lengths and deep energies. 

B. Calculation with non-universal corrections 

We first attempt to reproduce the previous results 
by including non-universal corrections to the Efimov 
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Figure 3: Ratio between the effective range r e and scatter- 
ing length a of the scaled LM2M2 potential, as a function 
of the scaling coefficient A. The vertical solid line indi- 
cates the physical value (A = 1). The dashed vertical line 
indicates the scaling where the scattering length diverges 
(f/a — > 0). The dashed curve shows the result based on 



the analytical formula ( 16 ) using the values of a and a as 
a function of A. 




Figure 4: Helium 4 dimer energy E as a function of scat- 
tering length a. To clarify the figure, these quantities are 
normalised to E e = h 2 jmr\ and r e , and raised to the 
power 1/4 and —1/2, respectively. The solid red curve 
corresponds to the dimer of the LM2M2 potential scaled 
by a coefficient A. The physical scattering length of helium 
4 (A = 1) is indicated by the vertical gray line. The dotted 
line represents the universal limit of small binding energy 
and large scattering length [Eq. ( 14 1 with a(E) — > a]. The 
dashed curve represents the two solutions of the effective 
range approximation of Eq. 1 14 1. The dotted-dashed curve 



corresponds to the results of the separable Gaussian po- 
tential (frj> . 



theory. For each value of A, we can determine the 
energy dependence a(E) of the scattering length of 
the scaled LM2M2 potential XV (r). The energy E 2B 
of the two-body bound state is then readily obtained 
from the equation: 
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Figure 5: Efimov spectrum for helium 4: trimer energy E 
as a function of the scattering length a. To clarify the fig- 
ure, these quantities are normalised to E e — ti 2 /mr 2 and 
r e , and raised to the power 1/4 and —1/2, respectively. 
The zero of the energy axis corresponds to the three-body 
threshold. The red solid curves are the energy curves ob- 
tained by scaling the LM2M2 potential. Both trimers and 
dimer (rightmost curve) are displayed. The vertical gray 
line indicates the scattering length corresponding to the 
unsealed potential, i.e. physical helium 4. The dashed 
curves correspond to the Efimov trimer spectrum accord- 
ing to the universal theory. The 3-body parameter is ad- 
justed to match the first excited trimer state. The straight 
dashed line corresponds to the dimer energy in the univer- 
sal limit of large scattering length. 



The low-energy expansion of a(E) up to the effective- 
range term (right-hand side of Eq. |5])) can already 
reproduce the realistic two-body energy for a > 2.5r e - 
see Fig. [4] However at this order of expansion there 
are actually two solutions to Eq. (14 1, the lowest- 



energy solution being unphysical. The two solutions 
merge and disappear at a = 2r e . The presence of this 
extra dimer is an artefact which completely changes 
the 3-body physics. Note that this problem does not 
occur for negative effective ranges, as in the case of 
narrow Feshbach resonances |41| . To remedy this 
problem, we consider an improved analytical expres- 
sion for a(E) that is obtained from the separable po- 
tential Eq. ([7]) adjusted to reproduce the scattering 
length and effective range of the scaled LM2M2 po- 
tential. Then there is only one solution to Eq. (14 1 
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and its energy matches very well that of the LM2M2 
potential - see Fig. |1J Interestingly, the agreement is 
even slightly better than the effective-range approxi- 
mation itself. 

Having set a(E) to properly describe two-body 
physics, we perform trimer calculations using Eq. @ 
with a fixed cutoff P whose value is adjusted to repro- 
duce the second 3-body dissociation point - the point 
where the first excited trimer reaches the threshold. 
One can see in Fig. [6] that the corrections bring some 
improvement, but the agreement with the LM2M2 
results is only partial. We then determined the re- 
quired variation of the cutoff P in order to obtain per- 
fect agreement with the LM2M2 results for both the 
ground and first excited trimer states. The variation 
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Figure 6: Similar plot to Fig. [5] Here again, the helium 4 
curves are indicated by solid red curves, and results taking 
into account two-body corrections to the universal theory 
are indicated by dashed curves. 



Figure 8: Similar plot to Fi g. |5| for the calculations with 
the gaussian potential Eq. (|6j). The dotted curves cor- 
respond to two-body interactions only, while the dashed 
curves correspond to calculation with an additional three- 
body interaction. 
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Figure 7: Three-body cutoff parameter P adjusted to get 
agreement with the LM2M2 results, as a function of en- 
ergy, for the ground-state trimer (dots) and the excited- 
state trimer (squares). One can see that their variations 
are inconsistent. 



is represented in Fig. [7] and is inconsistent for both 
states. To remove the inconsistency, we can assume 
more generally that P depends on both the energy 
and the scattering length, but no clear pattern arises 
from such considerations. It should also be noted that 
the required variation of P is dependent on the choice 
of the high-energy behaviour of a(E), in other words 
it is model-dependent. Thus, while it can be a practi- 
cal way to characterize non-universal observations, as 
done in Refs. |15 | Hl) l l2"2"]. it does not seem to be very 
meaningful. 



C. Calculations with a separable potential 

We then attempt to reproduce the 4 He results with 
a simple two-body potential V having the same low- 
energy properties as the LM2M2 potential. For both 
the Gaussian potential ([6| and separable potential ^ , 
we proceed as follows. 
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Figure 9: Similar plot to Fig. [5] for the calculations with 
the separable gaussian potential Eq. Q. The dotted 
curves correspond to two-body interactions only, while the 
dashed curves correspond to calculation with an additional 
three-body interaction. The symbols show the measured 
dissociation point of the ground-state trimer for different 
species. 



For each A, we adjust the parameters of the poten- 
tial so that the scattering length and effective range 
coincide with those of the scaled LM2M2 potential. 
As mentioned before, the binding energy of the two- 
body bound state then matches very well that of the 
LM2M2 bound state over a wide range of energies - 
see Fig. |4j 

We then calculate the three-body bound states with 
the adjusted potential. Remarkably, we already find 
relatively good agreement with the LM2M2 three- 
body states - see Figs. [8] and [9] The values differ es- 
sentially by an energy shift, which is not unexpected 
since the short-range three-body phase may not be 
properly set by these simple two-body interactions. 

We finally add a three-body interaction (also chosen 
to be Gaussian or separable), and adjust its strength 
in order to reproduce the second three-body dissocia- 
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tion point. We then find very good overall agreement 
with the LM2M2 calculations, as shown in Figs. [8] 
and El 



D. Role of the effective range 

The previous results clearly indicate that the non- 
universal variations of the trimer energies can be re- 
produced to a great extent from the knowledge of the 
effective range. The role of the effective range was 
pointed out before in a previous studies fl8 | fT9l . 130 1 132], 
In particular, Ref. [35] also looked at deviations from 
universality with finite-range potentials. This study 
considered the scattering length ajj gg for which the 
trimer dissociates. The relative variation of this quan- 
tity with respect to its value a^ gs in the universal 
spectrum was found to be: 



a diss a diss 
a diss 



= C 



V a Jdi 



diss 



(15) 



with C = 1.3 ± 0.4. However, we could not com- 
pletely confirm this formula. While we found the value 
C w 0.99 with the Gaussian potential, which is consis- 
tent with the results presented in Ref. [55] for the same 
Gaussian potential, the present calculation with the 
LM2M2 potential gives C « 0.58. The value there- 
fore seems to be somewhat dependent on the type of 
potential. 

It is interesting to note that the value of the disso- 
ciation point a^jgg itself, rather than its non-universal 
variation, was found experimentally to be in a narrow 
range. Experimentalists have measured the ground- 
state trimer dissociation point near several diver- 
gences of the scattering lengths [17], and for differ- 
ent species [H |3J [5] E0 EH [H HZ] , and found that it 
almost always occurs around a^iss ~ —9.9a, where 
a — 27r/r(l/4) 2 / v( j-yy is the average scattering length 
of van der Waals potentials [35] . Since there is an ap- 
proximate relation between a and the effective range 

tv nana: 



= 2 T(l/4) 4 a 
a 3 (2tt) 2 a 



a 

- - 1 
a 



(16) 



this corresponds to o,A; ss /r e ~ —2.76. Experimen- 
tal measurements of o,^ ss /r e for different species are 
represented in Fig. [9] It is quite remarkable that the 
scaled helium potential also gives a dissociation point 
a diss/ re ~ — 2.82 (or equivalently a^iss/^ = —10.26) 
which lies in the same narrow range. Furthermore, the 
model potentials considered in this paper (Gaussian 
and separable potential), without any 3-body inter- 
action, also give a consistent value a^ ss /r e ~ —2.70 



of the dissociation point. This shows that the use of 
these model potentials to interpret experiments leads 
to a good estimate of the dissociation point, as noted 
before [211 |2"2"] . Why it does so is however a puz- 
zling question, since we know from calculations with 
other or deeper potentials, that the dissociation point 
can change significantly [23]. This suggests that, un- 
der some condition yet to be discovered, the effective 
range may not only determine the non-universal vari- 
ations, but also the three-body parameter, thereby 
determining the full three-body spectrum. 



IV. CONCLUSION 



In this work, we have studied the Efimov physics of 
three helium-4 atoms, as a simple but realistic exam- 
ple to understand the non-universal variations of the 
trimer energy with respect to the scattering length. 
We found that non-universal two-body corrections to 
the universal theory alone are not sufficient to fully 
reproduce these variations, and a variable three-body 
parameter is needed. However the variations of this 
three-body parameter do not seem to follow any sim- 
ple rule, and are model-dependent. On the other 
hand, ad hoc but simple two-body potentials such 
as separable Gaussian potentials adjusted to have the 
same scattering length and effective range reproduce 
remarkably well the trimer energy non-universal vari- 
ations, and can constitute relatively accurate substi- 
tutes for single-channel realistic interactions provided 
that a small and localised three-body force is intro- 
duced to properly shift the energy. This indicates that 
the non-universal deviations which were observed in 
|151 fT5] and could not be explained by these effective 
potential models are likely to be due to more subtle 
effects such as multichannel coupling, non-trivial ef- 
fects of three-body forces, unless they simply result 
from underestimated uncertainties either in the scat- 
tering lengths or experimental data. 

Provided that such subtle effects are absent, our 
results also suggest that in general, beyond the usual 
Efimov universal scenario occurring for higher excited 
trimers, the whole spectrum follows a more specific 
class of universality determined only by the scattering 
length, the effective range, and the strength of a three- 
body localised force setting the three-body parameter, 
i.e. the position of highly-excited trimers. 

We thank F. Ferlaino, R. Grimm, P.S. Julienne, 
J. D'Incao, C.H. Greene and E. Braaten for helpful 
discussions. Some of the numerical calculations were 
performed on the HITACHI SR11000 at KEK. 
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